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ABSTRACT 



The QCD finite energy sum rules method is used to show that the parameter of the K^-K'^ 
mixing 13 is mainly determined by the value of gms{dG ^uf^^yd) and the vacuum expectation 
values of four-quark operators. Assuming the hypothesis of vacuum dominance and/or unitarity 
symmetry to estimate the latter, it is found that 5 = 1.2 ± 0.1 



^Copy of the paper published in Phys. Lett. B174 (1986) 104. Received 8 July 1985; revised 
manuscript received 11 April 1986. A couple of minor misprints is corrected. Comments, 
related to a (relatively small) change of the result due to changes in the values of incorporated 
phenomenological parameters, are added. 



The Kl-Ks mass difference is rather sensitive to the mixing of the t-quark and u- and 
c-quarks, so the analysis of the K^-K^ system provides us with useful information about the 
values of mixing angles in the Kobayashi-Maskawa model and about the phase parameter 5 of 
CP- violation. The K^-Ks mass difference Am can be presented as a sum of the long-distance 
dispersive contributions Am^ and the short- distance contributions Am^^: 

Am = Am^ + Am'^^ . (1) 

Am^^ is related to the matrix element of the effective AS = 2 hamiltonian (see e.g. Ref. 
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2Re(K"|/7f^}='|ir") 

{Gl/16n')Fix„e,)iM^/mK){K^\d\K''), 



(2) 



where O = {sLladLYois{^)~'^^'^ is the renormalization-invariant operator of the hamiltonian 
which arises in the calculation of the well-known box diagram [0, Gp is the Fermi constant 
and My/ is the mass of the l^-boson. The function F{xj, 6j) has the following form [Q]: 

Re[\lS{xc)r]i + \lS{xt)r]2 + 2\c\tS{xc,Xt)^^] . 



F{x,,ei, 



(3) 



Here x, 



are defined as 

S{x) 

S(^Xi, Xj) 



m'j/M^r, Xi = V*iVis {Vij is the Kobayashi-Maskawa matrix), and the functions S 
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The coefficients rji make allowance for the strong-interaction corrections in the leading logarith- 
mic approximation. For = 100 MeV they have the following numerical values: t]i = 0.7, 
ri2 = 0.6, and rj^ = 0.4 [Q. The coefficients Aj are connected with the mixing angles in the 
following way: 

ReAc ~ sjcj, ReXt ~ s^c^Sj; ReAcAt ~ slcjsl (5) 



where s,- 



sin 6j and Cj 



cos 6j, and the experimentally acceptable hypothesis that sin ^3^0 



and sin 02 ^ sin 9i is used. The mixing angles have the following numerical values: si = 
0.229 ± 0.006, ci = 0.9735 ± 0.0015, and 0.016 < S2 < 0.095 §. 

As follows from eq. (2), in order to calculate the short-distance contribution to the K^-Ks 
mass difference, it is necessary to find the value of the matrix element 



which is usually expressed through the dimensionless parameter B. 



(6) 



As for the long-distance contributions, they have been estimated using different phenomeno- 
logical approaches. For instance, the authors of refs. [^, |^ propose to estimate these contribu- 
tions by inserting the low-lying states between two AS* = 1 weak non-leptonic hamiltonians. 
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More precise estimates have been obtained in refs. P-||TO|; the corresponding values of the 
parameter D = Am^/Am are D = 0.10 ± 0.41 D = 0.33 ± 0.37 |[, and D = 0.46 ± 0.13 



A lot of attention has been paid recently to estimating the parameter B. The first such 
estimate -5 = 1- has been obtained in ref. |^ using the vacuum dominance approximation. 
Substituting this value into formula (2) and taking into account only the contribution of the 
c-quark with the mass = 1.3 — 1.5 GeV, we find that the short- distance contribution is about 
40-45 % of the total value of Am. The t-quark contribution is small: for rrit = 40 GeV and 
tb > 10~^^s it does not exceed 2% of Am. Thus at i? = 1 the short-distance contribution does 
not saturate the experimental Kl-Ks mass difference, which indicates that the long-distance 
effects may be also important. Unfortunately, the current estimates of B depend considerably 
on both the particular model used in the calculation and the values of various parameters 



involved in the model p|,[|ri[]-||l6l. Thus at present it is not clear what is the real value of this 
important parameter. 

This paper presents a new estimation of the value of B within the QCD finite energy sum 
rules (FESR) method, which has demonstrated its efficiency in studying the properties of low- 
lying hadronic resonances [|17[ |l^ . 



Within this approach, the problem to be solved is reduced to computing the following 
three-point correlator [ p!9[| : 

r,.(p, = .7 i.iv.Mw^ - nv) X mMlm« (7) 

Pfj.qyT{p^^ [p — qY, q^) + other structures 

at small q and large (euclidean) p: > 1 GeV^. Here = d^y^'j^s is the interpolating field 
of the K'^ meson 

{0\jlmK\p)) = tfKP, , /k = 1.17/. = 156 MeV 

The matrix element (6) is related to the amplitude T{t) = T{—t, —t, 0) by means of the disper- 
sion relation in the variable t: 

T{t) = J ds ^^'^^ — subtractions (8) 



, (K°|0|ir°) A 



where the single-pole contribution corresponds to the transition form factor of the meson, 
and the dots stand for the contributions of higher states. 



In the vacuum dominance approximation T^^,y assumes the form 

h 
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^i^a{p) = i / dxexp{ipx){T{f (x)sL{0)-fadL{0)))o 
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and the resulting value of B proves to be B^^ = 1. Thus, there remains to be computed only 
the function A^,^ = T^,^ — T^jf, which is responsible for all the departures from the vacuum 
dominance prediction for B. In other words, within the sum rule approach we will calculate 
the value B — 1. 

In the kinematical region we are interested in, there are three distinct contributions to T^^, 
viz: the perturbation theory (PT) contribution; the power corrections due to the non-zero 
vacuum expectation values (VEV's) of local operators pO]; and finally, the power corrections 



proportional to two-point correlators depending on the momentum q (or, in other words, the 
VEV's of bilocal operators |^). These bilocal power corrections are due to the fact that one of 
the external momenta (g) is small, and thus, the contribution from the region of large y ~ 
in the integral (7) cannot, generally speaking, be disregarded. 

It is easy to check that the leading PT contribution to T^i, is to be completely assigned to 
T^°. Thus, neglecting all the non-leading PT corrections we find that A^y|p2_,_oo = A^^^ + A^,^, 
with A^i^ (Aj^^) standing for the terms proportional to the VEV's of (bi)-local operators. 

The account of local operators with dimension < 6 leads to the following result for the 
function A^^: 

Aj;^ = p^" q" {-5{pq) (asG^) / 1927t'^ -A{dsds) (9) 
-A{ddss) + 2{ssss) + 2{dddd) 
+ms{gdG^u(^^ud) /2A'k'^)p~'^{p - q)'^ 
+other structures 

where the designation {dssd) stands for {di^'yaShShladh) and so on; = G'J^^^G'Jj^^,] G^u = GJ^j^t"; 
tr{tH') = i^'^^ and a^, = ^^[7^,7.]- 

To estimate the contributions due to bilocal operators one needs to construct the Wilson 
expansion for the T-product i J dx exp{ipx)T{j^{x)j^{0)) at large euclidean p (or, equivalently, 
at small x). A straightforward calculation shows that in the leading order in as, the bilocal 
contribution to the tensor structure p^gj, is suppressed by the factor p~^ and, thus, can be 
neglected within our approximation. 

To proceed, the combination of four-quark operators appearing in (9) proves to transform 
as a component of an SU-^{3) 27-plet. This means that the corresponding VEV is at least 
of second order in the unitary-symmetry breaking parameter. Moreover, the VEV of each 
four-quark operator in (9) vanishes if the vacuum saturation procedure is to work. Thus, 



the corresponding contribution to A^^^ is "doubly" forbidden and can be safely neglected. On 
the other hand, a straightforward estimation shows that should both of these suppressions 
be absent in the next-to-leading approximation, the corresponding contribution to AB might 
be uncomfortably large despite the loop suppression factor ~ Os/vr. To clarify the situation 
we have calculated the a<j-corrections to the coefficient functions of the four-quark terms in 
expansion (9) and have found that the appearing extra terms are die out after applying the 
vacuum saturation. 

In order to extract information on the value of B we employ the FESR technique to the func- 
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tion T(t) multiplied by (m|- + 1) to nullify the effect of the (unknown) single-pole contribution 
to the RHS of (8). The final sum rule has the form 



B - 1 



2 



(10) 



nis {gda^^G^yd)/lQTx'^m\ff 



K 



where 



(27ri)~iim 



T(— s — ze, —s — ie, 0) — T(— s + ie, —s + ie, 0) 



and So = 1.2 GeV is the duality interval of the K meson [18 



Up to now we have neglected all the effects due to the non-zero anomalous dimensions of 
the operators under consideration. In the leading log approximation the account of these effects 
is carried out without any difficulty. We define the renormalization-group invariant quantity 
B = i?(yu)[as(yu)]^^/^ and make use of the renormalization group technique p2| to rewrite 
eq. (10) as 

B = (1 + ms{gda^^G^^d)/l6n^m'^j^ff^\so)[as{so)]~'^^^ . 



To estimate the numerical value of B, let us reduce the matrix element (gdG^^cx^^d) accord- 
ing to the relation {gdG^yij^jjd)= ml{dd), = 0.8 ± 0.4 GeV^ I^S] and use PCAC relation 



(m„ + md){dd) 



^ r2 2 



whence 

B-l = -{f^ml/QATT^mlf^)2msml/{m^ + m,). 

The quark mass ratio is known with a relatively high accuracy [Q: 2ms/ (^u+^d) = 25.0±2.5. 
Finally, choosing A = 100 MeV we find 



B = 1.2 ±0.1 



(11) 



Prediction (11) is in agreement with the calculation of B within the lattice approach pSl , 



and also with the constraints from above B < 2.0 ± 0.5 [|I4| and the result B = 0.9-1.2 



Iq] , which were obtained by the method employing the dispersion representation for two-point 
Green functions. At the same time, our result is about three times as large as the value of 
B ^ 0.3-0.4 computed in the bag model |1^ and by using PCAC hypothesis and the S'f/'^(3) 
symmetry in ref. [|T3[] . 

To conclude, we have considered the problem of calculating the matrix element (6) within 
QCD FESR approach. We have found that the bulk of contributions to the relevant Green 
function T^^, keeps within vacuum dominance approximation, and thus, their effect on the 
value of B can be taken into account even without any calculations! In our opinion, this is 
the principal advantage of our method which most of the other approaches apparently do not 
have. It has allowed us to compute the parameter B with a fairly high accuracy in spite of 
the noticeable uncertainty involved in the determination of the non-factorizable contributions 
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to T^jy. Indeed, we have shown that in the leading order in a^, the deviation of the actual 
value of B from its vacuum dominance estimate B^^ = 1 is within the following limits: 0.1 < 
[13 — 1) < 0.3, where the main uncertainty is due to our poor knowledge of the matrix element 
{gdG^ya^yd). However, the overall magnitude of i? — 1 proves to be small, which gives the 
above margins for B. 

The authors are grateful to V.A. Matveev and A.N. Tavkhelidze for their constant help 
and support. We thank N.V. Krasnikov, Z. Maki, V.A. Rubakov and M.E. Shaposhnikov for 
numerous discussions. One of us (K.G.Ch.) is grateful to P. Langacker for a useful discussion 
of the problem. 

Note added. After this work has been issued as a preprint, we have learnt of another calcula- 
tion oi B with the result B = (0.33±0.09)[as(/i^)]^''^, which is quite different from ours. In fact, 
the authors of ref. have combined the effective chiral Lagrangian approach with the finite en- 



ergy sum rules method and calculated the parameter B{t = 4m|-) = (0|O|fC°-ft'°)/(2/|^m|^/3) 
rather than the parameter B defined as in the formula (6). 

To our mind, one of the possible reasons for this disagreement could be the noticeable 
variation of the function B{t) between t = Amj^ and t = (within the chiral perturbation 
theory the quantity 4mj^ can hardly be considered small in any way). It should be stressed 
that our approach meets no such problem since from the very beginning we are dealing with 
B{0) = B. 

Another possible explanation to the discrepancy could be the bad convergance of the power 



correction series for the correlator examined in ref. ||2^. Let us recall that in this work the 
account of the two first power corrections have led to the decrease of the result by five times, 
while in our analysis the power corrections to the (known beforehand) factorizable contribution 
does not exceed 20% of the latter. 

Leaving aside these technical subtleties, we feel that the main advantage of our approach is 
the exact account of the factorizable contributions and that the only chance to vary the result 
(11), say, by a factor of 2 or 3 is to have a strong violation of the vacuum saturation hypothesis 
for four-quark operators. 

At present the problem of estimating the accuracy of the latter hypothesis and the possibility 
of the exact account of factorizable terms within our approach are under consideration. 
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Comments 



Since the time the paper was pubhshed (1986) the need in an accurate numerical value 
for the Bx has considerably raced up. Precise experimental data on K'^ — system have 
appeared that allowed to reliably determine e'/e. The program of computing the next-to- 
leading correction to weak Hamiltonian in the effective theory approach has been completed 
[|l|. However, no significant improvement in calculation of the matrix element has been achieved. 
The reason is clear: it is a non-PT problem. As such it is now actively discussed within lattice 
approach but the corresponding results are still far from being perfect. 

Our calculation was made within the standard continuous QCD and the QCD Sum Rule 
Approach. The method we used factors out completely (and thereby takes into account) the 
factorizable contribition to the matrix element in all orders of PT and OPE expansion. The 
remaining non-factorizable part happens to be small in comparison to factorizable term, the 
latter being accounted exactly in our set-up. 

The method is open to improvement: both high order PT corrections and next order power 
corrections in addition to those computed in the paper can be added. Note that PT corrections 
are really mandatory to compute to match the available NLO coefficient functions (see, e.g. 

i). 

Having in mind a solid theoretical basis of our calculation we have just updated the pre- 

i f=3) 

diction for Bk by taking into account a significantly changed value for A^ from 100 MeV 
to the value of over 400 MeV, which corresponds to the world average value as{Mz) ~ 0.118. 
The updated version of eq. (11) with as{l.2 GeV^) = 0.69, estimated in the leading order, is : 

B = 1.0 ±0.1 

Note that the bulk of the change of the parameter B is due to the normalization (the change 
in the factor 0^(1.2 GeV^)^^/^) while other parameters of the calculations entering eq. (10) did 
not change much during last 15 years. 
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